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We demonstrate that entanglement of two macroscopic nanoelectromechanical resonators - cou- 
pled to each other via a common detector, a tunnel junction - can be generated by running a current 
through the device. This can be most efficiently achieved if the two oscillators are initially both 
prepared in their ground states. We propose two kinds of setups where the generation of entangle- 
ment can be realized by two different means. In the first setup, the oscillators are indirectly coupled 
via common fermionic reservoirs with long coherence times. While this setup gives valuable insight 
in the physics of this open quantum system, the second proposed setup, an Andreev entangler, 
represents a novel and feasible way of entangling two nanomechanical oscillators. In the Andreev 
entangler, a split Cooper-pair that coherently tunnels to each oscillator, mediates their coupling and 
thereby generates entanglement between them. 

PACS numbers: 03.65.Ud, 85.85.+J 



Introduction. — Since the early days of quantum me- 
chanics, entanglement has been one of its most fascinat- 
ing key elements pQ. The study of entanglement not 
only is concerned with fundamental questions but also 
with applications, for instance, entangled qubits being 
the building blocks for quantum computers. Two of the 
main questions when studying entanglement are: How 
can one generate entanglement in a multipartite system 
and what is a good measure for entangled states? For 
discrete variable multipartite systems with, for instance, 
qubits as representatives, the former two questions have 
been extensively studied theoretically as well as experi- 
mentally [2H3- 

Only within the last few years, nanoelectromechani- 
cal systems (NEMS) have been providing an experimen- 
tal and theoretical playground for studying quantum ef- 
fects on a macroscopic length scale. Among the latest 
experimental achievements are cooling a nanomechani- 
cal resonator to its ground state [5HZ] and engineering of 
non-trivial quantum mechanical states [8]. On the theo- 
retical side, entanglement generation in NEMS has been 
a goal for many years. In the pioneering work of Ref. [3], 
a route towards entanglement in NEMS was proposed 
making use of a global non-adiabatic change of the inter- 
action strength in a one-dimensional chain of nanome- 
chanical oscillators, leading to entanglement of two spa- 
tially well separated oscillators of the chain. This type 
of entanglement has, however, up to date not been re- 
alized in experiments. Subsequent proposals in (hybrid) 
nanooptomechanical systems have suggested to generate 
entanglement between a moveable mirror of a cavity and 
a collective spin ensemble of an atomic medium [10 , a 
Bose- Einstein condensate jTT], [T2] , the cavity field itself 
[T5] . or another mirror [T4TU6| . 

Entanglement creation and verification of quantum 
systems, as for instance, atoms and photons is a theoret- 
ically well understood subject and routinely done exper- 
imentally in various labs [I7HI5] . Entanglement of quan- 



tum systems with constituents being more than just one 
atom or photon, as for instance, collective spin ensembles 
and Josephson qubits, have also been studied recently 
[2T)1 |2"T] . These kinds of systems are often times discrete 
variable systems with a finite dimensional Hilbert space 
and true quantum systems in the sense that they do not 
have any classical counterpart. With the laws of quan- 
tum mechanics describing atomic and subatomic systems 
and the laws of classical Newtonian mechanics describ- 
ing macroscopic objects, there must exist a boundary be- 
tween the quantum and the classical world. This quan- 
tum to classical boundary is still subject in contemporary 
physics [221 [53]. Emerging questions are: What causes 
the quantum to classical transition? Where is the limit 
that causes objects to behave quantum mechanically and 
not classically? How macroscopic can an object be to be 
still considered quantum? 

This work is motivated by exactly these questions. Our 
aim is to investigate the quantum mechanical feature per 
se, entanglement, in a quantum harmonic oscillator that 
behaves as classically as it gets (in the sense of the Ehren- 
fest theorem). We like to study the challenge to gener- 
ate entanglement between two nanomechanical oscilla- 
tors, being representatives of two mesoscopic continuous 
variable systems. To this end, we propose two funda- 
mentally new schemes towards entanglement generation 
in NEMS. 

The first scheme, Setup (A), consists of two nanoelec- 
tromechanical oscillators coupled to a common measure- 
ment device/reservoir and is depicted in Fig. [lja). This 
setup being rather academic serves us as a model and 
provides us with new insight into the situation where 
two nanomechanical oscillators are effectively coupled 
through two common fermionic baths. This is a less 
studied case compared to the more established case of 
coupling two oscillators to a common bosonic bath, for 
instance, within a Caldeira-Leggett model [24]. For a 
more feasible version of Setup (A), we introduce a second 





FIG. 1. (Color online) (a) Schematic setup of a bipartite continuous variable quantum system, realized as two nanoelectrome- 
chanical oscillators (green) in a tunnel junction setup which share common fermionic reservoirs (blue). Yellow lines indicate 
tunnel junctions that are needed for the generation of entanglement, (b) Two nanomechanical oscillators are effectively coupled 
in an Andreev entangler setup due to a Cooper-pair that is split. Each of the two electrons of the Cooper-pair is assumed to 
tunnel from the superconductor (orange) onto a different oscillator. 



system, Setup (B), where we bypass the common (right) 
reservoir. Setup (B) is inspired by an Andreev entangler 
setup [25] and schematically depicted in Fig.Qb). 

Method. — Both of our proposed setups are best an- 
alyzed in a system-bath approach. The total system is 
described by the generic Hamiltonian 



H = H S + H B + H. 



SB 



(1) 



where H$ and Hb describe the system and bath degrees 
of freedom, respectively. The system Hamiltonian for 
Setup (A) and (B) is given by H s = j=i,2 Pi / 2m * + 
rriiQfx^/2 which corresponds to two nanomechanical os- 
cillators with effective mass mi, frequency Qi, and po- 
sition and momentum operators Xi and pi, respectively. 
For simplicity, we later assume two identical oscillators, 
i.e., fli — D,2 — and m\ = 7712 = m. The cou- 
pling between the system and the bath is given by H$b 
which we will specify below for Setup (A) and (B) sep- 
arately. Our goal is to investigate the time evolution of 
entanglement between two nanomechanical oscillators. A 
Gaussian bipartite system is fully characterized by its 
covariance matrix CTy(t) = TVs[ps(i){£j, £y}/2], where 
Ps(t) — Tr B[p(t)] is the system's density matrix and 
£ = (xi,Pi,X2,p2) T ■ As a measure for entanglement, 
we use the logarithmic negativity Em [17]. The loga- 
rithmic negativity is directly related to a(t) — cr Tl (t) 
through its smallest symplectic eigenvalue A_ (t) and given 
by Ejtf(t) — max{0, — log 2 [2A_ (t)]}, where a{t) is the par- 
tially transposed (with respect to subsystem 1) covari- 
ance matrix [23]. We calculate the time dependence of 
the entries of a(t) with an equation of motion for the 
system's density matrix ps(t), where we treat the non- 
Markovian regime by employing a time convolutionless 
master equation method |29| . The non-Markovian dy- 
namics implies that the system's time evolution depends 
on its history, reflected in the equation of motion by 
time dependent damping and decoherence kernels. In 



the Schrodingcr picture, the master equation in the Born 
approximation reads (we put K = 1) 

ps(t)=-i[H s ,ps(t)} (2) 

- / drTvB [H SB ,[H SB (T-t),p s (t)® p B }} ■ 
Jo 

With Eq. ([2]), expectation values of arbitrary system op- 
erators X (or products of them) can be calculated by 
solving d/dt(X(t)) = Tr s [ps(t)X]. 

Setup A. — The first setup we propose, can be seen as 
a natural extension of the experimental realization of one 
oscillator coupled to a single atomic point contact [30] . 
Setup (A) consists of two nanoelectromechanical oscilla- 
tors that are both coupled via an atomic point contact 
to common electronic reservoirs, cf. Fig. [ija). Impor- 
tantly, the oscillators are not directly coupled to each 
other (as e.g. in Ref. [26]), but only through the elec- 
tron reservoirs which act as fermionic baths. The yellow 
tunnel junctions in Fig. [lja) are sufficient to generate 
an entangled state of the two oscillators. The Hamil- 
tonian for Setup (A) is of the form of Eq. (fll), where 

Hb = J2 r £ r t Pr' l l ; r + J2i £ i W ' ^ 'i refers to the two leads 
(depicted blue in Fig. HTa)). The system-bath coupling 
is meditated via the tunneling Hamiltonian 



HsB = ^2 ("f0 + Ix&nWllPr + H - c - 



(3) 



l,r n=l,2 



where for small oscillation amplitudes we approximate 
the tunneling matrix elements as linear in the oscillator 
displacement, and for the sake of simplicity assume all 
tunneling amplitudes to be real and take the coupling to 
be symmetric. In order to treat non-Markovian effects 
of the fermionic reservoirs properly, we need to include 
an energy dependent density of states in the leads. This 
is done with a Lorentzian shaped density of states which 
means that an electron with energy Ei in the left lead can 
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tunnel into states of the right lead with energy e r , broad- 
ened by L c [24] I3TH33] . If we compared Setup (A) with a 
scenario where the right reservoir is cut into two halves 
and therefore each oscillator is coupled to independent 
reservoirs on the right hand side, there would be no pos- 
sibility to have any coherent effective coupling between 
the two oscillators. Hence, the common right reservoir is 
crucial for the generation of entanglement in this setup. 
The equation of motion for the reduced density matrix 
of the system up to second order in the tunneling can be 
written as (details are given in Ref. [34 ) 



T=0, V=0.1,y,=0.1 



Ps(t) 



H s + l JC^ ) (t)(x 1 +x 2 ) 2 lPs (t) 
(t) [xi + x 2 , [xi + x 2 ,ps{t)]] 

IC+ (t) [&i + X 2 , [pi +P2,ps(t)]] 
iC^\t) [ Xl + X 2 ,{Pl+P2,Ps(t)}} 



0.04 
tq 0.03 
0.02 
0.01 




(4) 



FIG. 2. (Color online) Logarithmic negativity for L c = 0.3 Q 
(red) and L c = 0.6 Q, (blue). Here, the temperature T = 0, 
V = 0.1 f2, and 7^ = O.lVmQ. As can be seen, the qualitative 
dependence on the parameter L c is rather small. 



We now investigate the system's dynamics, where we 
take the vacuum state as an initial state, which is Gaus- 
sian, allowing us to use the logarithmic negativity as a 
measure of entanglement. Here, we only discuss the case 
of zero electronic temperature. We show in Fig.[2]the log- 
arithmic negativity for j x — and its dependence 
on L c . We clearly see that the initially separable state 
becomes entangled right after the interactions have been 
suddenly switched on. Due to this quench, energy conser- 
vation can be violated for a short time, which is the cause 
for the entanglement. Experimentally, the switching in 
this tunneling setup can be achieved by gates control- 
ling the tunnel coupling via a resonance. Physically this 
means, that the rise time of the gates has to be shorter 
than l/f2, the timescale of the oscillators. Typical rise 
times of electronic gates can be as short as 60 ps [35] . For 
an oscillator with a resonance frequency of ~ 500 MHz 
(which can be brought into its ground state at low tem- 
peratures) these rise times are sufficient to accomplish 
the sudden switching. In its subsequent evolution, the 
entanglement oscillates in time, before it slowly decays. 
We find that low bias voltages show an increased loga- 
rithmic negativity compared to high bias voltages (not 
shown in the figure). 

Setup B. — The major drawback of Setup (A) is the 
assumption that the which-path information of the tun- 
neling electrons is lost because the common right reser- 
voir cannot spatially distinguish the electrons. To cir- 
cumvent this assumption, we now introduce an alterna- 
tive realization based on a superconducting reservoir in 
the center and two independent normal metal reservoirs 
which nonetheless is capable of generating entanglement. 
This realization, Setup (B), is of an Andreev entangler 
type. In its original sense, the Andreev entangler [25] 
relies on the splitting of a Cooper-pair which is a spin- 
singlet. There, the split Cooper-pair from an s-wave 
superconductor can coherently tunnel via two different 



quantum dots into different leads. Throughout this pro- 
cess, the spin-singlet is preserved. Therefore, the An- 
dreev entangler is capable of generating nonlocal spin- 
entangled electrons. In our context of NEMS, we use 
only charge properties of the split Cooper-pair as a me- 
diator for an effective coupling between the oscillators. 
We show in Fig. |T|(b) what such an Andreev entangler 
setup could look like. A superconducting island (or- 
ange) serves as a source of Cooper-pairs which can tunnel 
onto two different (conducting) nanomechanical oscilla- 
tors. The process where the Cooper-pair is split and 
one electron tunnels to the lower oscillator and the other 
one tunnels to the upper oscillator gives rise to an effec- 
tive coupling between them. The major improvements 
of Setup (B) over Setup (A) are: (i) the right common 
reservoir of Setup (A) is now cut into two halves; the re- 
sulting two normal metal leads are entirely independent 
(ii) both normal metal leads see the same tunneling pro- 
cesses; therefore their independence does not result in 
which-path information (Hi) phase coherence lengths of 
superconductors are typically much longer than of nor- 
mal metals which allows for a larger spatial extend of the 
center reservoir. Processes where one Cooper-pair tun- 
nels to the upper/lower lead are not taken into account 
since they won't lead to any entanglement, but rather to a 
background tunneling current. As for instance described 
in Ref. [25], such processes are also energetically not fa- 
vorable due to interaction effects in the leads and / or tun- 
neling via a dot with a finite charging energy. Another 
possibility to favor the mentioned tunneling process, is to 
couple the superconductor to two Luttinger liquid leads 
as described in Ref. [36] . In NEMS this could be achieved 
by modeling the oscillators as one-dimensional nanowires 
or carbon nanotubes. We focus here on transport at en- 
ergies in the superconducting gap A, which is the largest 
energy scale of the problem. 
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The system is again described by the Hamiltonian in 
Eq. ([IJ. Now, the bath consists of three independent 
reservoirs, two normal metal leads, and a central su- 
perconducting region. The bath Hamiltonian becomes 
H B = Hj + H V + H D with H a = E^c^LAfc. 
(q = U,D) and = Eto^L^- Here E k = 
\J (ek — Ps) 2 + A 2 is the quasiparticle energy and ftko- 
is the quasiparticle annihilation operator in the super- 
conductor. We take the superconductor to be grounded 
and each lead to be held at bias voltage V. The system- 
bath interaction is mediated by the tunneling of a split 
Cooper-pair. The process we focus on is therefore de- 
scribed by the effective tunneling Hamiltonian 



H. 



k± .k2,ks,k4 



^U^D^D-k 1 i^Uk 2 -t C I-k 3 i C Ik 4 ^ 



(5) 



+T u T D b Dki1 .b u _ k2 ^c I _ k3 ^c Ik ^ + H.c. , 

where T a = "foa+lxa x a an d the electron annihilation op- 
erators for spin f and l are given by c^-y — Uk/3kf+v k /3_ k ^ 

and c-ki — u kP-kl — v k p^, respectively. The effective 
Hamiltonian, cf . Eq. ^ , exhibits two important but dif- 
ferent contributions. The first contribution stems from 
terms of order 7o7x- These terms are qualitatively sim- 
ilar to terms of order "f x in Setup (A), meaning that to 
second order in j x we expect a qualitatively similar be- 
havior for the logarithmic negativity. Only considering 
these terms in the master equation, all kernels IC(t) and 
JC{t) for Setup (B) can be calculated in the same way 
as the kernels for Setup (A) (for details on the equa- 
tion of motion for ps(t) for Setup (B) see Ref. [34]). In 
Fig. |3j we present the results for Ejsr(t) for Setup (B). 
From this figure, we already see the similarity between 
Setup (A) and (B) and conclude that it is possible to 
entangle two nanomcchanical resonators which are cou- 
pled to three independent reservoirs within an Andreev 
cntangler setup. The second contribution of the Hamil- 

T=0, V=0.1 , 70=03 , y x =0 .05 




to a direct coupling between the two oscillators, and are 
a new element in the Andreev entangler setup. Interest- 
ingly, these terms can even generate entanglement in the 
Markovian regime. In the low bias limit, we can apply the 
rotating wave approximation since the oscillators cannot 
be excited by the applied bias and excitations can only be 
swapped between them. The resulting master equation 
is of Lindblad form and given by 



Psit) 



i[H s ,ps{t)] 



(6) 



Ap s (t)A- -{AA, Ps (t)} 



where the Lindblad operator A = 



it; 



+ a-^al, and 



p = PuPdPsPs with p a and ps being constant den- 
sity of states in lead a and the superconductor, re- 
spectively. We also introduced bosonic annihilation 
flj = (xiy/mfl + ipi I V rafi) j 'v2 and creation operators 

it 



a\ = (xi^mVL — ipi/ \/mf2)/v2- To exemplify the dy- 
namics generated by the Lindblad operator A, we choose 
|ni = \.ri2 = 1) as an initial state (n% being the occu- 
pation number of oscillator i). Due to Eq. (|6| the dy- 
namics will only generate states |ni = 2, ri2 = 0) and 
| rii = 0,ri2 = 2). The system will stay in this three- 
dimensional Hilbert space. The initial state evolves to 
pi{t) |1, 1) +pz(t) (|2,0) + |0,2)). A definite indicator for 
entanglement in this case is the prefactor P2(t) in front 
of the Bell state. Here, we just want to make the point 
that within the Andreev entangler setup we can generate 
an entangled state even in the Markovian regime. 

Conclusion. — We introduced two setups, where we 
have shown that it is possible to generate entanglement 
between two spatially separated nanoelectromechanical 
oscillators, even if these two oscillators are not directly 
coupled to each other. In Setup (A), the two oscillators 
are indirectly coupled via two common fermionic baths 
which should have a rather long coherence time. With 
Setup (B), we introduced an entirely new way of entan- 
gling two nanomechanical oscillators in an electric setup 
based on the working principle of an Andreev entangler. 
Here, the coherent process where a Cooper-pair is split 
leads to an effective coupling of the two nanomechanical 
oscillators. For both setups, the generated entanglement 
persists over many oscillator periods, before the system 
becomes separable again. 
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FIG. 3. (Color online) Logarithmic negativity for L c — 0.3 
(red) and L c = 0.6 (blue). Here, the temperature T = 0, 
V = 0.1 fi, 70 = 0.3, and j x = 0.05Vmfi. 

tonian Eq. (15J) is due to terms of order 7^ which give rise 
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I. DAMPING AND DECOHERENCE KERNELS OF SETUP (A) 



Here, we give details on the equation of motion for ps(t) and the appearing kernels for Setup (A). We first simplify 
the system-bath coupling Hamiltonian and write 



H SB = S " E + 5 ™# f » 



(S.l) 



n=l,2 



where, S n and E denote arbitrary operators acting only on the Hilbert space of system and bath, respectively. We 
define the system and bath operators as S n = 70 + lxX n and E — ipjipn respectively. With this, the equation of 
motion for the reduced density matrix of the system up to second order in the tunneling, Eq. (2) of the main text, 
can be written as Eq. (4) of the main text, where we defined the time dependent damping and decoherence kernels 



4 x) w = 


[ dr 
Jo 


K (1 \t) 


+ K<»( 


*?>(*) = 


[ dr 
Jo 


K {1 \t) 


+ K^( 




[ dr 
Jo 


'k^\t) 




ic ( ^\t) = 


[ dr 
Jo 




_#(2)( 



sin(f2r) , 



sin(f2r) 



with 



K^(t) = ^{{E(t),E\o)}}, 
^ (2) W = ^([^)^ f (o)]>. 



(S.2) 
(S.3) 



We see from the equation of motion for ps(t) that its time evolution is governed by the time dependent damping 
and diffusion kernels K, l ^J_{t) and K+?(t). We will briefly sketch the calculation of the kernels X 1 / 2 (i). Due to the 
fact that the equation of motion for ps (t) is of second order in Hs b > the only possibility of including non-Markovian 
effects is by considering an energy dependent density of states in the leads. Including non-Markovian effects leads 
to a finite correlation time in the leads, and is a key ingredient for the entanglement. Non-Markovian effects in the 
reservoirs and their influence on entanglement of two quantum systems have also been studied in Refs. [IrE]. The 
kernels K<~ m \t) are given by 

Jf("»)(i) = i J d £l J de r J(e l7 e r ) e^'-^[n( £l ){l - n(e r )} - (-l) m n(e r )ll - n( £l )}] , 

where n(e x ) — [e 13 ^*^^^ + is the Fermi distribution function and we introduce an energy dependent spectral 
function 

J{si,e r ) = ^25{ei - e k )S(e r - e q ) . 
k . (1 

To account for a finite lifetime of quasiparticles in the leads, the (5-functions are smeared out and replaced by 
Lorentzians of width L r 



J{e h e r ) 



E 

fe . q 



(ei-e k ) 2 + L 2 c {e r -E q y + Ll 



(S.4) 
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The largest contribution of each of the independent sums in Eq. (S.4) will come from the energies close to the Fermi 



level of each lead. We further restrict ourselves to the regime of low applied bias voltages (V < L c ). This is just a 
simplification in order to keep the number of parameters in the problem as low as possible. With these assumptions, 
the energy dependent spectral function can be written as 



1 



(ej - S r ) 2 + LI 



(S.5) 



Physically, Eq. (S.5 1 means that an electron with energy ei in the left le ad c an tunnel into states of the right lead with 



1 /2 

energy e r , broadened by L c [H[SHZ]. With the spectral function in Eq. (S.5), the time dependent kernels K,f,_{t) and 



'1/2 

£/,_(£) can now be calculated analytically, but the resulting expressions are too lengthy to be stated here. 



II. EQUATION OF MOTION AND TIME DEPENDENT KERNELS OF SETUP (B) 

The equation of motion and the calculation of the kernels for the master equation in the case of Setup (B) are very 
similar to Setup (A). Therefore we will only briefly sketch the main ideas. 

The quasiparticle operators of the Hamiltonian Hj are related to electron annihilation operators through the 
Bogoliubov transformation 

Cfet = u k /3 k f + v k pl H , 
c-ki = - Vfe/^t ' 

where Uk — \/l/2 + ^fc/(2£ , / c ) and Vk = \/l/2 — £fc/(2_Efc) with = Sk — ps- As before, we introduce system and 
bath operators 

50 = 7oc/7on , 

51 = lODlxUXu , 

52 = iQUlxD^D , 

53 = IxUlxDXDXU , 

and 



E — ^ D D-k 1 \J D Uk 2 \ C I-k 3 i c Ik i \ ) 



(S.6) 



respectively, which allow us to write the system-bath coupling Hamiltonian in a similar way as in Eq. (S.l). With 
this, the equation of motion for Setup (B) is similar to Eq. (4) of the main text and can be written as 



Ps{t) 



-i[#s + i/C (2) (i)(£i+x 2 ) 2 , (0 s(i) 
- K>l\t) [xi + x 2 , [xi + x 2 ,ps(t)]] 

+ /C< f 1) (<) [Xl+X 2 , [Pi +p2,p S (t)]] 

+ iC^\t) [x 1 + x 2 ,{p 1 +p 2 , Ps (t)}} 



K 



:3,i)/ 



+ (*) 
~ic ( + 2 \t) 



xi£ 2l [x\X2, ps(t)]] 
xix 2 , \p1P2, Ps(t)]} 
x 1 x 2 , [xif> 2 +x 2 pi,p s {t)]] 
±ix 2l {xix 2 ,ps(t}}] 

XlX 2 ,{piP2,Ps(t)}] 

xix 2 , {xip 2 + x 2 pi,ps(t)}] 



(S.7) 
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where we have defined the time dependent memory kernels in Eq. (S.7) as 

rt 



/d 2) (t) 
/ci 2) (t) 



/C<_ 3 ' 2) (i) 



/d 3 ' 2) (i) 



/C^' 3) (t) 



dr 



(i 



= / 



= / dr 



dr 



(It 



= / dr 



= / dr 



dr 



(It 



= / dr 



K (2 



)(r 
)(r 
)(r 
)(r 
)(r 
)(r 
)(r 
)(r 
)(r 
)(r 



The kernels i^ 1 / 2 )^) are given in Eq. (S.2) and (S.3 



-K<? 
-K<- 2 



(-r 

(-T 
(-T 
("T 
(-T 
(-T 
(-T 
(-T 
(-T 



7o7* cos(fir) , 



2 2 

7o7 2 

.,2. ,2 



sin(fir) . 



%lt cos(Or) , 



sin(57r) . 



7o 2 7 2 



r y x cos(rir) cos(fir) . 



7 



m 2 £l 2 

4 



sin(fir) sin(fir) . 



cos(rir) sin(f2r) 



7^ cos(fir) cos(Or) . 



7 



TO 2 f2 2 

4 



sin(fir) sin(17r) 



7 



77 ? 



cos(rir) sin(rir) 



with i? given in Eq. (S.6). The actual calculation of the kernels 



for the Andreev entangler setup goes along the same lines as for setup (A). 
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